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The procedure of the gauge fixing is an essential part of QCD |L| and 
however final results do not depend on the gauge, different forms of gauge 
conditions are useful in different settings of physical problems. For example, 
in high energy scattering in QCD the axial gauge has proved to be useful [Q], 
while in the OPE analysis the Fock-Schwinger |3J (sometimes called the 
coordinate or radial) gauge was applied (for discussions and derivation see 
g] and also @). 

In another physical situation, where the time axis is singled out, as e.g. 
in the heavy quarkonium theory, the modified coordinate gauge can be 
convenient. This gauge was used recently in the context of equations for the 
quark || and gluon || Green's functions, displaying the property of chiral 
symmetry breaking and confinement. 

There is another set of studies where an emphasis is made on formulation 
of gauge theory without gauge-dependent degrees of freedom from the very 
beginning, and the role of dynamical variables is played by 2-forms [|TI] 



loop variables [|Tl] or some auxilary (vector - type) variables [|12j. These 
completely gauge invariant approaches encountered their own difficulties and 
as a matter of fact many gauge-invariant observables are easier to calculate 
using gauge dependent diagrammatic rules. 

Both in the coordinate gauge J| and in its modified form j7j the shape of 
the contour C(x), in the integral, connecting vector potential and the field 
strength, 

A p (x) = J dz u a Pfl (z)F up (z) (1) 

C(x) 

is fixed and consists of straight lines. Inessential for physical results, it may 
be inconvenient in the course of computations. In particular, in the confining 
phase of QCD, when the QCD string is formed between two colour charges it 
would be advantageous to choose the contours C lying on the world sheet of 
the string; in this case one could do simplifying approximations as in |8|, [J, 
namely to keep only Gaussian field correlator. The decoupling of ghosts, 



known to occur for the gauges ([!]) (see [O], and references therein) is 
also an attractive feature, which suggests to look for generalizations of (Jl|) 
with arbitrary contours C. 

There is also a deeper reason for the interest to the gauges of the type 
([ID. It lies in the fact, that (0) can imply dynamical connection between field 
variables and geometrical (i.e. contour) variables. We briefly outline this 
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possibility at the end of the present paper. 

The gauge condition of the type we are interested in was introduced 
for the first time in [|TJ[]. In the present paper we give a refined treatment 
of this gauge, paying special attention to some important details, missing 
in the original paper. Let us briefly mention them. To define this gauge 
condition correctly, the set of contours C, determining the gauge must satisfy 
some additional requirement (eq.(|3|) of the present paper). This condition is 
essential for the representation (|l|) to hold true. With this requirement we 



are also able to formulate the gauge condition in the local form (eq. (|TTD of the 
present paper). An immediate use of the generalized contour gauge which 
also has not yet been discussed in the literature is the ability to give a short 
and direct proof of the nonabelian Stokes theorem [^, |TB[ as we do below 
in this paper. Another line of development pursued below is the derivation 
of the QCD Lagrangian in terms of the dual vector potential and contour 
variables. 

Let us proceed with the definition of the generalized contour gauge. Let 
M be a (i-dimensional connected Euclidean manifold. We choose some sub- 
space Mo, Mo C M which in general may be disconnected and of lower 
dimension than d. In the simplest case, to be considered below, M consists 
of the only one point Xq. 

For each point x G M \ M we define the unique smooth contour , 
Xq G M connecting points x and Xq. The contours are parametrized as 
follows: 

CI : = z»(s,x); s G [0, 1]; z»(0,x) = x "; 2f(l,x) = x» (2) 

The map M \ M — > M defined above is naturally extended to M — > M 
by setting C%° to be the unit contour: z^(s,x ) = Xq^. The resulting map 
M — > M is assumed to be smooth. In the particular case when M consists 
of the only point xq it means, that the manifold M should be contractible. 
This requirement plays an essential role in what follows. 

Let us choose two arbitrary points z^s.x) and z^s'.x) on the given 
contour C in such a way that the point z^{s',x) lies between points z^(s,x) 
and Zfj,(l,x) = x^ (if s is natural parameter, it simply means that s < s'). 
We assume the following condition - for any s, s' there exists s" such that 

z^(s,x) = z,i(s",z(s',x)) (3) 
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The geometrical meaning of (|3]) is simple: for any point z lying on some 
contour C* its own contour C| Q coincides with the corresponding part of the 
contour . The eq.(^) does not mean, generally speaking, that contours 
and from different points X\ 7^ x<i have no common points except Xq. 
The condition that contours C* should not selfmtersect (the only condition 
discussed in jTJ|) is necessary but not sufficient to guarantee (0) and therefore 
to derive (§) and (|TT|) below. The defined set of contours forms an oriented 
tree graph without closed cycles according to @. 

Let us now start with the gauge potential A p (x) taken in some arbitrary 
gauge and perform the gauge rotation 

Afa) = n + (x)A ll (x)n(x) + - n + (x)d IM n(x) (4) 

where 



J-0 



and integration goes along the contour C% . The important point is the 
differentiation of the phase factors [|K| which is a well defined procedure 



for our choice of contours since the function z p (s } x) is given. The contour 
derivative reads: 

d^(x) = niA.umr) - i 9 n(x)A p (x Q ) dr[){J ' ] 



dx p 



dz u (s, x 



+ ig J ds Zu ^ X a Pfl (z)U(x,z(s))F pu (z(s))U(z(s),x ) (5) 



where 



/ \ dzp(s,x) . . 

a P»( z ) = ( 6 ) 

By Xq(x) in the second term at the r.h.s. of fl5|) we denote the initial point 

Xq for the contour with the end point x. This term is trivially absent if M 
consists of the only one point. 

Substituting (|) into the (£|) one gets: 

A',(x) = A p (x ) ^+ 
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+ U M J d s 



dz„(s,x) 
ds 



a pft (z)U(x,z(s))F up (z(s))U(z(s),x ) (7) 



Taking into account the condition (|3]) and the gauge transformation prop- 
erty U(x Q , x)U(x, z(s))F pu (z(s))U(z(s), x ) — > F' (z(s)) we arrive to the final 
result 



A'(x) = A p (x ) 



8xq(x) 
dx„ 



1 

+ J ds 



dz v {s,x) dz p (s,x) 



ds 



dx, 



KM S )) (8) 



This formula was proposed in |14j] and used without derivation in ||. In 
the rest of the paper we take xq to be a unique point for all contours and 
therefore dx^/dxp = 0. 

The eq.(f|) leads to important local condition for vector-potential. To this 
end, note, that solving (|3|) with respect to s' we find: 

s' = f(s,s" } x); f(s,s,x) = 1 

Substituting s' = f(s, s", x) into (|3|) and differentiating with respect to s one 
gets: 

dz^s, x) _ dzp(s", z(s', x)) dz p (s', x) df(s, s", x) ^ 



ds 



ds' 



ds 



dz p (s',x) 

Then putting s' to be equal to unity and multiplying both sides of (H) by 
t M (x) = (dz fM (s,x)/ds) s=1 we get 



ds 



dz u (s,x) dz p (s,x) (dzp,{s,x) 



ds 



dx,. 



ds 



KM*)) 



8=1 



ds dzJ^X± dz^x) ! = Q 



ds 



ds 



(10) 



where g(s, x) = (df(s, s", x)/ds) s ii =s . The second equality holds by virtue of 
(JD and the third due to antisymmetry of F' . 
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The condition (|T0| ) can be easily understood taking into account, that 
phase factors along the contours C% specifying the gauge are equal to unity: 



Xo 



Since ([II]) holds for all x, one gets: 

A'^x) t,(x) = (11) 

Specific examples of the gauge condition discussed in the present article 
are known in the literature. These are the radial or Fock-Schwinger gauge 
H and its modified forms |7|, §]. Note that due to the topological restrictions 
stated above these gauge conditions are always defined in some neighbour- 
hood of their origins but might not be well defined globally, in particular, in 
the case of topologically nontrivial M. This was noticed in different respect 
also in . 



As an illustrative example let us consider the use of the generalized gauge 
condition for the nonabelian Stokes theorem. There are different proofs of 
this theorem in the literature [T5|, ITBJ, but what we are going to present is 



perhaps the simplest one. It is close in spirit to the paper ||15|| . Namely, we 
define the gauge condition in such a way that potential A^(x) on the contour 
is expressed as a function of field strength F^ u {u) defined on the (arbitrary) 
surface, bound by the contour. Then rewriting gauge- invariant Wilson loop 
in this gauge we obtain a relation, valid in the chosen gauge and as the last 



step put it into gauge-covariant form. It was done in [15] for the completely 
fixed axial gauge condition, which is a convenient choice in two dimensions 
(or for planar surfaces in higher dimensional case). Our procedure allows one 
to choose an arbitrary surface S bound by the simple contour C = dS and 
therefore the gauge condition we use entirely depends on the shape of S. We 
parametrize the surface by w^(s, t); s, t 6 [0, 1]. and choose an arbitrary point 
Xq^ on the surface in such a way, that u;^(0, t) = x^. If s = 1 then u> M (l, t) 
goes along the contour C and w^(l, 0) = w^(l, 1) according to dC = 0. 

The following important remark is in order. It is usually assumed that S 
has the disk topology, and the contour C = dS is unknotted, in this simplest 
case we are free in our choice of Mo, which may consist of only one point, 
what we actually have used. For this topology it is always possible to define 
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a set of contours obeying (|^) by, for example, continuous deformation of the 
planar disk with the radial contours. But in topologically nontrivial cases the 
proof should be modified, (see in particular ]r7|). In the rest of the present 



paper we are concentrating on the disk topology. 

According to the gauge potential A^(z) is related to F fiu (z) in the 
following way: 

AM°A) = hs' dZ ^f )) d -^^KMs',t)) (12) 
o ' s XfM ^ ' 

Equation ([12]) is actually nothing else than the Stokes theorem in its infinites- 
imal form. It is well known that the generalization to finite contours is non- 
trivial in the nonabelian case, in particular the integral J s F^da^ depends 
on the surface even if the contour C = dS is closed. But this integral does 
not enter by itself in the nonabelian Stokes theorem. Instead the quantity 

which should be considered here is a P-ordered exponent Pexp(ig J A^dx^ 1 ). 

c 

Substitution of fll~2|) into the definition of the P-exponent leads to the 
expression: 

ig J A l _ t (x)dx l _ l 

Pe O 

oo 

= 1 + £ (ig) n / •• / da, u (w (1) (s 1 ,t 1 ))..da p ^ n \s n ,t n )) 

71=1 J J 

F p(P (w (n) (s n ,t n ))..F^(w (1 \s 1 ,t 1 )) 9 > t 2 > ..>t n ) (13) 

Note the ordering procedure in (|T3|) - only the points along the contour C 
are ordered with respect one to another, i.e. ordered in parameters t iy while 
the integrals over Sj are taken independently for each ti. 

To bring ([l3j) to the gauge covariant form we introduce phase factors along 
the s-direction on the surface, which are equal to unity due to (|TT|), i.e. we re- 
place F^(w(s,t)) -> G^(w(s,t)) = U(x ;w(s,t))F^(w(s,t))U(w(s,t);x ) 
If the point x does not lie on the contour C the gauge-covariant answer 
reads: 

ig J A„dz„ ig J da^G^Az) 

Pe w =U(x*,x )Pe s U{x ,x*) (14) 



where the meaning of the ordering simbol V is explained in (0). Under the 
gauge rotations both sides of (0) are transformed in the same way. The 
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more often used gauge-invariant form of ([TJ|) is 

ig J A^dzp. ig J da fiv (z)G llv (z) 

TrPe c =TrVe s (15) 

We stress again, that the exact meaning of the symbol V is completely de- 
termined by the choice of the set of contours, defining the gauge which may 
be done as the most convenient one for a given application of the nonabelian 
Stokes theorem. 

It is worth noting that all formulas used above hold true in abelian case 
too, the ordering operations are not necessary in this case. In particular, 
for the given abelian gauge field strength F^ v (x) one can define a set of 
contours {C(x,xo)} obeying all necessary conditions discussed above and 
obtain a set of the gauge potentials {A^x)} by using formula (||). Then all 
field configurations {A^(x)} are gauge equivalent as it is clear from (£|). In 
other words one can say, that variation of the contour's shape in (|8|) with the 
fixed in abelian theory leads only to gauge variation of the field A^{x) 
and therefore the contour degree of freedom is not dynamical in this case, 
all dynamics is encoded in the field strength tensor. The situation changes 
however if we introduce monopoles into the theory. The relation between 
field strength tensor and gauge potential takes the form: 

F^(x) = d^A u (x) - duA^x) + 2ngJ2G IIU (x) (16) 

where 

G>„(x) = e^afs / da aP (y)5 {4l \x - y) 
s 

and £ is an (arbitrary) surface, which is a worldsheet of the Dirac string. 
Then the contour C{x) defining the gauge rotation Q(x) in @ can be smoothly 
deformed until it intersects E. If such intersection happens, gauge potential 
A^{x) receives a singular contribution proportional to 

J dz v J da afS (y) e upaf s5 {4) (z - y) 

It is seen, that even though Dirac strings are gauge degrees of freedom, 
the contours C may become dynamical ones in a theory with monopoles. 
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We leave this set of questions for the future publication and now focus our 
attention on the dual form of the QCD partition function. 

There are different ways to integrate out gauge potentials in nonabelian 
theories (see for example |12|], fl5l , |T8| and references therein). The gauge 



condition discussed in the present paper looks promising for that purpose 
because it gives an explicit expression for in terms of F^ u . It should be 
noticed, that (||) indeed singles out a unique representative for any gauge 
orbit. To see this, let us consider two gauge-equivalent connections A and 
A = uj + Auo + i/g uo + duj and perform rotations @)\ 

A 1 = Q + [A]A Q[A] + -fl + [A]dfl[A] 

A' = n + [A] AQ[A] + ^-n + [A}dQ[A] 

where integration in Q[A] and Q[A] goes along one and the same set of 
contours. One easily finds, that A'{x) = uj+(xo)A'(x)uj(xq) where xq is the 
base point, i.e. resulting potentials coincide up to the global gauge rotation. 
Let us proceed with the complete gauge fixing. Our strategy is a straight- 



forward extension of We choose the 4-bein p\(x)\ fM,i — 1..4 smoothly 
depending on x: 

A set of orthogonal curves x^ = x^s 1 ), i = 1..4 is defined in such a way that 



The potential A^x) = A % {x) p^(x) is taken to vanish on one of the chosen 
lines; assuming this line to be defined by the conditions x % = const, % = 1..3 
one gets: 

A^xis 4 )) = A^x 1 ,x 2 ,x 3 ) = 
The piecewise contour z^(s,x) is then defined as follows: 

3 

Zfl (s,x) = xl + Y: fds%(x(s*)) (17) 
=i 



"0 
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The path goes along the lines x l consequtively with respect to i, i.e. initial 
points in any integral from (|T7| ) coincide with the end points for the previous 
integral: = x(sq). It leads to the following local conditions for the 

potential: 

A 3 = A\x 2 ) = A 2 (x\ x 3 ) = A 4 (x\ x 2 , x 3 ) = (18) 

The eq.(^|) with the contour's definition ( |T7|) allow to express potential in 
terms of the field strengths. 

In the simplest particular case when p 1 = 8 l , the expression (|18|) defines 



the so called completely fixed axial gauge, used in [|15|]. One can easily check 
that the following conditions on the usual coordinate components of the 
vector potential hold true: 

A 3 = Ai(z ) = A 2 (xq,z ) = A^x ,y ,z ) = (19) 

It is seen that it is a particular case of generalized contour gauges discussed 
in this paper. 

For completely fixed gauges the following relation is very useful 



J DA5(CGF)5(G - F[A}) = 8{I[G)) (20) 

where I[G] is the Bianchi form I[G] = D^G^ V and G is the dual tensor 
G^u = ^/^^-^uapGap. Delta function 5{CGF) fixes gauge completely. The ex- 
pression (p0|) is local and Lorentz-covariant, taking into account that locally 
flTjp is equivalent to (|I9~D we conclude, that (|20| ) works for the generalized 
completely fixed gauges as well as for axial gauge. It is straightforward then 
to integrate out vector potentials A^ in the partition function. The QCD 
partition function reads: 

Z[J] = J DADipD4>5(CGF) 

/DG5(G-F[A])e /d4 K-^ G - +J - G - + ^^ [A]+m) ^) (21) 

where the field strength tensor -F M „L4] = ~ d u A^ — i[A^A u ], covariante 
derivative -D M [^4] = <9 M — iA^ and J^ v - external tensor source. 
It is convenient to introduce a nonlocal operator 

K^(x,y) = j ds^j^± d -^f^8^{z-y) (22) 
o S Xfl 
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Then the gauge fixing S(CGF) implies due to (^) that 

A,[G}(x) = J d 4 yK^(x,y)G up (y) (23) 
By using ( p0|) the partition function reads: 

Z[J] = J D^DkJ DG jJd^kUd^+f^A^G^) 

J ^x^-^G^+J^G^+^iiy^+m^+^^A^G}^ 

After the Gaussian G-integration one has: 

Z[J\ = J Dk^Di, (detA)-* el ' ^mi^+m)i,) 

where 

Kp{x) = ^ud^K{x) + J d 4 w^(w)i^t a il;(w)Kf(w,x) (26) 
and the propagator 

Af ppa (x, y) = -Lda^pJix -y) + if abc e^k c 7 (x)K^(x, y) (27) 

Using abelian language the second term in the r.h.s. of ([27]) represents the 
"monopole current", while k p plays the role of the dual vector-potential. 

In the absence of quarks and external currents the "dualized" partition 
function reads: 

Z = [ Dk^e-^ffW (28) 



The effective dual action S e ff[k] is similar to the one obtained in [15 



Seff[k] = 1 -TrlnA + J d'x J d 4 y [k] (x) ( A" 1 ) ^ p / b pa [k] (y) (29) 

and notice that dual field strength is a linear function of dual potential k in 
our formulation: ^ 

F a p[h] = -e aj 3^ v {d^k v - d v k^) 
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The stationary point k^x) of the effective action can be found from the 
solution of the equation 

5S eff [k) 



5k(z) 







(30) 



k=k 



In the large N c limit this yields in particular the effective action for quarks 
of the following form: 



S n 



d 4 u 



d 4 w U(u)i^t a ip(u) (kA' 1 



, ab 



)( 



_ (31) 

and the kernel KA 1 K is to be evaluated at the stationary point fc M . 

The dual action is essentially nonlocal, which is encoded in A(x, y). This 
property is common to dualized nonabelian theories (see |TB[). But in our 
case we have an additional freedom - to choose 4-bein p 1 ^ in such a way, that 
the operator K(x, y) takes the most convenient form for a given dynamical 
problem. This way of research which is ideologically close to the abelian 
projections method will be discussed elsewhere. 

Acknowledgments 

The work was supported by the RFFI grants 96-02-19184a, 96-15-96740, 
and partially by RFFI-DFG grant 96-02-00088G. The authors are grate- 
ful to A. Levin, M.Olshanetsky, V.Novikov for valuable discussions and to 
G.Korchemsky for pointing authors attention to the Ref.jbl . 



References 

[1] L.D. Faddeev, "Gauge Fields: An Intro to Quantum Theory", Addison- 
Wesley, 1991. 

F. Yndurain, "The Theory of Quark and Gluon Interactions", 2-nd ed., 
Springer- Verlag, New- York, 1993. 

[2] L.V. Gribov, E.V. Levin, M.G. Ryskin, Phys.Rep., 100 (1983) 1. 

[3] V. Novikov, M. Shifman, A. Vainshtein, V. Zakharov, Nucl.Phys.B249 
445 (1985). 



12 



[4] V.A. Fock, Sov.Phys.,12 (1937) 404. 

J. Schwinger, Phys.Rev.,82 (1952) 684, see also "Particles, Sources and 
Fields", vols.l and 2, Addison- Wesley, Reading, 1970 and 1973. 

M. Shifman, Nucl.Phys.B173 13 (1980). 

M.S. Dubovikov, A.V. Smilga, Nucl.Phys., B 185 (1981) 109. 

I.Balitskii, Nucl.Phys.B 254 (1985) 166. 

Yu.A. Simonov, |iep-ph/9704"30Tl 
Yu.A. Simonov, |hep-th/9712250] 



Yu.A. Simonov, |hep-ph/9712248| 

S.Mandelstam, Phys.Rev. bf 175 (1968) 1580; 
Ann.Phys. 19, 1 (1962) 25. 

Yu.M. Makeenko, A.A. Migdal, Nucl.Phys.B188 (1981) 269, 
for a review see A.A. Migdal, Phys.Rep.102 (1983) 199. 

P. E. Haagensen, K. Johnson, Nucl.Phys.B439 (1995) 597. 

S. Leupold, H. Weigert, |hep-th/96040T5 . 



S.V. Ivanov, G.P Korchemsky, Phys.Lett. B154 (1985) 197. 

S.V. Ivanov, G.P. Korchemsky, A.V. Radyushkin, Sov. J. Nucl.Phys. 44 

(1986) 145. 

M. Halpern, Phys.Rev.D19 (1979) 517. 

Y. Aref'eva, Theor.Math.Phys. 43 (1980) 353; 
N.Bralic, Phys.Rev. D22 (1980) 3090; 
Yu.A. Simonov, Sov. J.Nucl.Phys. 50 (1989) 134; 
M. Hirayama, S. Matsubara, |hep-th/ 97 12120] . 



M. Hirayama, M. Kanno, M. Ueno, H. Yamakoshi, frep-th/9806~098 



H. Reinhardt, |hep-th/9608T9T| . 

F. Lunev, J.Math.Phys. 37 (1996) 5351. 



13 



